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We consider strong interactions at finite density in mean field theory, through an effective la- 
grangian that can describe both nuclear matter and quark matter. This lagrangian has three cou- 
plings that are all fixed by experiment and no other parameters. With increasing baryon density 
we then find the following hierarchy. At nuclear density and above we have nuclear matter with 
chiral spontaneous symmetry breaking (SSB), followed by the pion condensed quark matter, again 
with chiral SSB, albeit with a different realization and finally a transition to the diquark CFL state 
which also has chiral SSB (and colour SSB), with yet another realization. To one's surprise at zero 
temperature (in mean field theory), at any finite density chiral symmetry is never restored! 

We find another remarkable feature and this is that the tree level mass of the sigma particle, 
that is set by experiment to about 800 MeV, has a crucial and unexpected influence on the physics. 
Strange quark matter and strange stars are ruled out for a sigma mass above 700 MeV and neutron 
stars with magnetic pion condensed cores, that could provide magnetic fields of neutron stars, exist 
only for a small interval, between 750-850 MeV, for the sigma mass. 



I. INTRODUCTION 

Neutron stars have been a subject of abiding inter- 
est for several decades, almost since Landau suggested 
their existence, shortly after the discovery of the neu- 
tron. There are a variety of astrophysical phenomena 
that arise from the physics of neutron stars. The su- 
pernova explosion through which the star is born, is a 
spectacular luminous event. Many neutron stars work 
as pulsars, which generate beamed radiation in their in- 
tense magnetic fields. Neutron stars in binary systems 
may accrete matter from their companions, giving rise 
to some of the brightest X-ray sources in the sky. Some 
neutron stars with super-strong magnetic fields produce 
occassional strong bursts of gamma rays (Soft Gamma 
Repeaters). Most of these phenomena require us to un- 
derstand the physics of matter at very high density, which 
govern the mass and the size of neutron stars. In other 
words, one needs to have a clear understanding of the 
equation of state of the ground state of superdense mat- 
ter. Although much effort has gone into this enterprise 
over the last four decades it still remains poorly under- 
stood. Why? 

Central densities of neutron stars are high, ~5 to 10 
times nuclear density p nuc — 0.17 fm~ 3 . For a single 
species, neutrons, this naively translates into a fermi gas 
with typical fermi momentum, fcy ~ 700 MeV. On the 
other hand nucleons have structure and a typical size of 
the order of a fermi (200 MeV) -1 . It is clear that at such 
high densities nucleons (neutrons) cannot be treated as 
elementary. They are composite and resolved. They are 
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colour singlet bound states of three valence quarks. At 
such high densities, therefore, treating nucleons as point 
particles interacting via two body (or more) forces will be 
inadequate. Yet most available equations of state adopt 
this approach and therefore fail to capture the correct 
physics. 

On the other hand, if we use quarks as the elementary 
degrees of freedom, we are presently bound by the fact 
that only perturbative calculations can be done for QCD. 
This implies that calculations can be done in QCD only 
at very high density when the theory is approximately in 
an Asymptotically Free (AF) phase. However, at inter- 
mediate and low density (close to nuclear density), where 
a nucleonic description is valid, we cannot use perturba- 
tive QCD as the coupling becomes strong and the physics 
non perturbative. This is the dilemma. 

There are attempts to model the physics by a 
two phase structure - a quark matter core with a 
hadronic/nucleonic exterior shell and crust. Since there 
is no simple way to link the two phases without using 
separate parameters for both, this description is some- 
what arbitrary. Further, the nature of the quark matter 
state is not clear - for example, is it in a spontaneous 
chiral symmetry broken state. 

Can we find a theory that can describe both these do- 
mains? We present, here, an Effective Chiral Intermedi- 
ate Lagrangian, L, that has quarks, gluons and a chiral 
multiplet of \n , a] that flavor-couples only to the quarks 

This yields the nucleon as quark soliton - a bound state 
of quarks in a solitonic background of scalar/pseudoscalar 
field expectation values (EV-s) that follow from the spon- 
taneous breaking of chiral symmetry 0, • I* 1 this way 
we can generate nucleon matter from these nucleons at 
low density with a transition to quark matter at high 
density but with the same effective Lagrangian covering 
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both domains. Such a unified description has not been 
given before and depends just on the coupling constants 
of the theory and not on any parameters. 

To begin with let us consider the two main features of 
the strong interactions at low energy. These are i) that 
quarks are confined as hadrons and ii) chiral symmetry 
is spontaneously broken (SSB) with the pion as an ap- 
proximate Goldstone boson. There is no specific reason 
that these two phenomena should occur at an identical 
temperature scale, though QCD lattice simulations show 
that for SU 2 (L) x SU 2 (R) they are close. The problem in 
giving an unequivocal answer to this question is that we 
are yet to find a solution to the non-perturbative aspects 
of QCD. 

Let us consider QCD with a two flavour SU 2 (L) x 
SU 2 (R) chiral symmetry. First let us address the ques- 
tion of what occurs at a lower energy scale, confinement 
or chiral symmetry breaking. The problem is that though 
there is a bonafide order parameter for chiral symmetry 
breaking - the mass of the constituent quark, the Wilson 
loop is no longer an order parameter for confinement, in 
the presence of dynamical quarks. However, by looking 
at the energy density or specific heat we can get a fair 
idea of the change in the number of operational degrees 
of freedom or particle modes. Such lattice calculations 
indicate that the change from the large number of degrees 
of freedom in the quark matter phase to few degrees of 



freedom in the hadronic one takes place in one broad step 
in temperature, indicating that the two transitions may 
be close for SU 2 {L) x SU 2 {R) ■ 

Also, if the chiral symmetry restoration (en- 
ergy/temperature) scale was lower than the confinement 
scale we would expect hadrons to show parity doubling 
below the confinement scale but above the chiral SSB 
scale. This is not seen in finite temperature lattice sim- 
ulations. 

Actually, QCD can have multiple scales 0- Apart 
from a confinement scale and a chiral symmetry restora- 
tion scale we also have a compositeness scale for the pion. 
The above considerations suggest these scales are respec- 
tively in ascending order in energy (temperature). 

Let us consider the interacting fermi liquid of nucle- 
ons - nuclear matter. As the baryon density is raised 
beyond overlap, we expect a transition to quark matter. 
An interesting question arises: Is the quark matter in a 
chiral SSB state with constituent quarks or is it, as is 
usually assumed, in a chirally restored state with current 
quarks? As we will see the Lagrangian L, given below, 
answers this question [(J Q ■ 

On the other hand there is some evidence for this in- 
termediate Chiral Lagrangian that has simultaneously no 
confinement but chiral SSB. Such an effective Lagrangian 
has quarks, gluons and a chiral multiplet of [jf, a] that 
flavor couples only to the quarks. 



L = -\(P^G% 9y{° + 1> - \{d,af l -(d^f ^ 2 (<r 2 + t? 2 ) - ^ (a 2 + vr 2 ) 2 + const (1) 
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The masses of the scalar (psuedoscalar) and fermions fol- 
low from the minimization of the potentials above. This 
minimization yields 

At 2 = -A 2 < a > 2 (2) 

It follows that 

to 2 = 2A 2 < er > 2 (3) 

Experimentally, < a > = f n , the pion decay constant. 
This theory is an extension of QCD by additionally 
coupling the quarks to a chiral multiplet, (if and a) 

alalia 

This Lagrangian has produced some interesting physics 
at the mean field level |j, |8| 

(i) It provides a quark soliton model for the nucleon 
in which the nucleon is realized as a soliton with quarks 
being bound in a skyrmion configuration for the chiral 
field expectation values (EV) 0, @ . 

(ii) Such a model gives a natural explanation for the 
'Proton spin puzzle'. This is because the quarks in the 
background fields are in a spin-isospin singlet state in 
which the quark spin operator averages to zero. On the 



collective quantization of this soliton to give states of 
good spin and isospin the quark spin operator acquires a 
small non zero contribution [jj. 

(iii) Such a Lagrangian also seems to naturally produce 
the Gottfried sum rule [Ic| . 

(iv) Such a nucleon can also yield from first principles 
(but with some drastic QCD evolution), structure func- 
tions for the nucleon which is close to the experimental 
ones |Tlj . 

(v) In a finite temperature field theory such an effec- 
tive Lagrangian also yields screening masses that match 
with those of a finite temperature QCD simulation with 
dynamical quarks 0] . This work also does not show any 
parity doubling for the hadronic states. 

(vi) This Lagrangian also gives a consistent equation 
of state for strongly interacting matter at all density 0, 

am. 

This L has a single dimensional parameter, f^, that 
is the pion decay constant, and three couplings, <?3, the 
QCD coupling, g y , the Yukawa coupling between quarks 
and mesons, that will be determined from the nucleon 
mass and the meson-meson coupling, A, which, for this 
model, can be determined from meson meson scattering 
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|l4j . No further phenomenological input will be used. As 
it stands, there is no confinement in this model, but may 
be dynamically generated as in QCD. 

Since this is posited as an effective Lagrangian, we 
should have an appoximate idea of its range of valid- 
ity. We find, somewhat in analogy with the top quark 
(large Yukawa coupling) composite higgs picture, that 
we can get a compositeness scale for the scalars in this 
model by using Renormalisation Group (RNG) evolution. 
We find that the wavefunction renormalisation for the 
scalars is inversely proportional to the running Yukawa 
coupling and thus naively vanishes when the Yukawa cou- 
pling blows up. For our theory such a ballpark scale falls 
between 700-800 MeV [13 ■ 

An independent and quite general approach in setting 
a limit to the range of validity of non asmptotically free 
(e.g. Yukawa) theories is the vacuum instability to small 
length scale fluctuations (or large momenta in quantum 
loop corrections) that plagues these theories, discovered 
and analytically proved by one of us 16]. The scale at 
which this occurs is of the same order as above. This is 
not very surprising since it is connected to non-AF char- 
acter of the Yukawa coupling [liil Il7| . This underscores 
the impossibility of doing loop calculations, unless we in- 
troduce a cut off, even though our L is renormalizable! 

Given these facts we shall use this theory at the Mean 
Field level to look at different phases of this field theory. 
To do this we must first establish the ground state of 
the Baryon Number B=l sector of this theory, i.e. the 
nucleon. 

The plan of the paper is as follows, In Sec. 2 we review 
the description of the nucleon in this model, which fixes 
g y . This is followed in Sec. 3 by a pedagogical description 

I 




The term V(a, n) = ^(a 2 + t? 2 - f 2 ) 2 is the potential 
functional. For the vacuum sector (no fermions) this is 
the quantity to be minimized. The minimum of V(a, 7?) 
occurs at a 2 + ir 2 = f 2 . This is the equation of a 3 sphere 
and thus allows for a continuous degeneracy. However, 
normally the choice, < a >= J^, < if >= is made, 
which is consistent with the pseudoscalar, < 7? >, having 
a zero Vacuum Expectation Value (VEV) to avoid spon- 
taneously violating parity in hadronic scattering. Also, 
with this choice the goldstone pseudoscalar excitations 
about the ground state are the right ones - the psue- 
doscalar pions. Once this choice for the vacuum state 
is made, it is clear that the vacuum state spontaneously 
violates the chiral symmetry, it changes under a 0(4) 
rotation, even though H does not. 



of how we may look upon nuclear matter at density above 
nucleon overlap. We then look at some phases of 2 flavour 
quark matter in Sec. 4 and point out that the phase with 
lowest ground state energy is the pion condensed phase 
with chiral SSB. In Sec. 5 these results are generalized 
to 3 flavours. In Sec. 6 we review and compare the pion 
condensed phase with the diquark condensed colour su- 
percondcting phase. We end this section with a review 
of the phase diagram of QCD at all density. In Sec. 7 
we remark on the phases to be used in constructing neu- 
tron stars in a following paper and discuss the validity of 
Mean Field Theory (MFT). 



II. THE NUCLEON IN THE CHIRAL LINEAR 
SIGMA MODEL WITH QUARKS 

The basic fields are the three component (isospin) pion 
fields, 7?(r), the scalar field, cr(r), which together form a 
real 4 component chiral multiplet that transforms as a 
representation of 0(4), the rotation group in four dimen- 
sions. The fermionic fields are the quarks which trans- 
form as a fundamental representation SUl(L) x SUr(R) 
which is isomorphic to 0(4). In this we ignore the gluon 
fields which are vector and therefore not expected to 
carry expectation values in MFT. However, corrections 
due to one-gluon exchange interactions between quarks 
can be included, but we shall not do so below, as their 
effect is small. We also neglect the pion mass corrections. 
The Hamiltonian which is invariant under the above 
group transformations then reads 



+ 0(m») (4) 

I 

We now move to the description of the ground state 
for the the baryon number B = 1 sector. 

1) The usual pattern of symmetry breaking is the space 
uniform VEV corresponding to the choice above. 

< cr >= /„., < 7? >= (5) 

This choice gives a spontaneous or Yukawa mass to the 
quark, m q = g < a >. 

The lowest energy for B = 1 sector, which corresponds 
to the quantum numbers of the nucleon, is M = 3m q = 
3^/tt - simply the mass of three quarks. 

2) Let us now consider the case when a fermion bound 
state can arise from a time independent expectation value 
(EV) that is locally space dependent. Clearly, to have fi- 
nite energy for this state requires that asymptotically (far 
away from the localized fermion source), the EV revert 
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back to the VEV above. We now move on to the so called 
skyrme hedgehog configuration, 

< a >= f n cos 9(r), < tt >= rfn sind(r) (6) 

where 6(r — > oo) = 0, from the finite energy condition 
and 9(r — > 0) = — tt, for the pion field to be well defined 
at the origin. 

It may be pointed out that we have chosen a 2 + tt 2 = 
f 2 . More generally, we need not fix this magnitude and 
vary a and tt independently. Given such a configuration, 
we solve the Dirac eigenvalue equation for the quark in 
the background field, 9{r). Due to the linking of space 
with internal isospace, neither the angular momentum, J, 
nor the isospin, /, commute with the Dirac Hamiltonian 
but only the sum K = J + I docs. K can then be used to 
label the eigenstates. The lowest, K = 0, valence state 
is a spin isopin singlet of the form 

iPk=o = Mr)\K = 0) (7) 

where, 

|X = 0> = -L(| + 1/2)|T) - |-1/2)U)) (8) 

the arrows designate spin and the halves the isospin. The 
eigenvalue equation is 

(-idicti + /3g y {<7 + i75*?.7f))^i<- = o = £o^/f=o (9) 

The \K = 0) state is a bound state, where R has the 
interpretation of the width of the potential. For a partic- 
ular profile, 9{r) = —ir(l — r/R), for r < R, and 6(r) = 0, 
for r > R, the eigenvalue dependence on R is sketched in 
Q(a),H. 

The nucleon may then be obtained as a bound state of 
three coloured quarks. Notice that the only degeneracy 
of the K = state is in colour, since this state is a 
spin isospin singlet and in making the nucleon we have 
exhausted this degeneracy, yielding a colour singlet. The 
energy of the nucleon state is given by 

E B=1 [e(r)} = 3e o [0(r)] + / d 3 x\m*f + (c^cr) 2 ] (10) 

in terms of a general profile, 6{r). This must be mini- 
mized with respect to variations of 0{r) to get the ground 
state energy, E mm . Finally, since we have eigenstates of 
K, this quark soliton must be projected into good states 
of spin and isospin. The mass, M, of this soliton with 
quark bound states which has the quantum numbers of 
the nucleon depends just on f n , g y and A. The depen- 
dence on A is marginal, so the only parameter that is free 
is the Yukawa coupling, g y . We fix, M = M nuc i eon . 

Actually, for such a nucleon, which is a colour singlet 
bound state of three valence quarks in a skyrme back- 
ground, in a linear sigma model, a generally accepted 
value for g y = 5.4 0(b). This is when the tt and a fields 
are varied independently, without making any ansatz, 



and the quark soliton so obtained is projected to give 
a nucleon with good spin and isospin. 

Let us now compare this to the energy of 3 free quarks 
in the uniform phase. We find that, 

M/(3g y U) < 1 

indicating that in our sigma model the nucleon is indeed 
a soliton with quark bound states. This analysis makes 
a further significant point. The threshold value of g y at 
which M I (3g y fTr) is first less than 1, is g y = 4. Therefore, 
for values of g y greater than 4, the quark soliton is always 
of lower energy than three free quarks and is thus bound. 
Besides the mass of this nucleon falls with increasing g y , 
as this controls the strength of the attractive potential in 
which the quarks bind. The conclusion, therefore, is that 
there is a maximum mass |l<Sl for the fermion in such a 
theory. 

Here we remark that we now have a determination for 
all parameters of our L: f v , g% and g y , leaving only one 
parameter, A to be set. As we shall show this can be de- 
termined from low energy meson-meson scattering data. 

Recently, Schechter et al made a fit to scalar chan- 
nel scattering data to see how it may be fitted with in- 
creasing i/i (centre of mass energy), using chiral per- 
turbation theory and several resonances. They further 
looked at this channel using just a linear sigma model. 
Their results indicate that for yfs < 800 MeV, a reason- 
able fit to the data can be made using the linear sigma 
model with a tree level sigma mass above but close to 
800 MeV. This sets the value of A. 

This completes the determination of all the parameters 
of our L, which is able to describe both nucleon and quark 
phases of dense matter. 

III. PHASES AT FINITE DENSITY 

We now move to the main theme of this work, which 
is to look at ground state of strongly interacting mat- 
ter at finite density. In our model all the phases con- 
sidered are characterised simply by different patterns of 
spontaneous symmetry breaking. As we shall see i) the 
nucleonic phase is characterised by the pattern of SSB 
in skyrmions ii) the Lee- Wick phase has space uniform 
SSB, with only the sigma field expectation value and iii) 
the pion condensed phase is characterised by a stationary 
wave, with a space dependent periodic variation in the 
sigma and neutral pion field expectation values. 

A. The nucleon or nuclear phase 

In this phase, as the name suggests, the quarks are to 
be found as bound states in nucleons. 

We found that the single isolated solitonic nucleon has 
lower energy than the three free valence quarks. Thus 
at low density we have a fermi gas of interacting nucle- 
ons. Long range potentials like pion exchange and tensor 
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exchange can be found by the appropriate two nucleon 
ansatz, in analogy to the two-skyrmion problem. These 
match well with the usual nuclear physics two-nucleon 
potentials It is then more reasonable to solve the 

nucleon many body problem by mapping on to the nu- 
clear many body calculation - since it is very difficult to 
solve a many body quark soliton problem! 

We know quite well the ground state of nuclear matter 
till, say, roughly two times nuclear density. From there on 
we can only model this phase at higher density, when the 
solitons start to overlap and their motion is obstructed, 
as a 'crystal lattice' of solitons 13j. We can then use 
the Wigner-Seitz approximation to convert the problem 
to a single cell problem. Scaling the size of the single cell 
scales the baryon density, since we have one nucleon per 
cell. 

It is appropriate to note that with such simplifying ap- 
proximations it is not fair to expect a realistic description 
but a pedagogically useful one. 

The Wigner Seitz approximation 

The crystal ansatz implies a picture where the solitons sit 
in a close packed hexagonal configuration. Each soliton 
is in a skyrmion chiral symmetry broken configuration 
for the [a, tt] fields that acts as a potential that supports 
a single bound state with a colour degeneracy of 3. The 
\K = 0) bound state is saturated for each soliton making 
the soliton into a colour singlet baryon. It is clear that 
when the solitons are close packed, the quark wavefunc- 
tion will leak out so as to minimize the energy. From the 
Kronig-Penny model it is expected that bound states will 
form a band with the number of states in the band being 
3iV, where N is the number of baryons stacked together. 
It is further known that the band will splay evenly about 
the single bound state energy, with the bottom of the 
band below and the top above the single soliton (K = 0) 
bound state energy. Due to the saturation of the bound 
state, the band will be completely occupied and a good 
approximation to the median energy of the band is the 
bound state energy of the single soliton above of radius 
R. Thus the energy per baryon can be approximated in 
the crystal configuration by the single soliton energy but 
with the parameter R setting the volume occupied by a 
single soliton. Clearly, the baryon density is set by R. 



At finite density the VEV becomes an EV and we can 
allow for 

< a > 2 + < 7T > 2 = F 2 (12) 

where F is now a variational parameter to be set by min- 
imizing the free energy. It is clear that this ansatz is 
rather restrictive. 

Note that at finite density there is just one constraint, 
namely 

i) < 7? > must vanish at the origin of the soliton for 
the pion field to be well defined. 



However, only for calculational simplicity we also main- 
tain ii) and iii), below: 

ii) F is space independent 

iii) < (j >= — F at r = 0, or 

< a >= F at r = R, or 

< 7? >= at r = R 

We choose a simple profile function, 0(r) — ir(l — -^), 
where R is the half-length of the cell. 

We may now calculate the Wigner-Seitz single soliton 
energy by calculating the quark energy eigenvalue eo> me- 
son gradient energy (second term in eq. (|13|l ) and sym- 
metry (condensation) energy (third term in eq. (|13|l 'l in 
this background: 

E b = 3e (6(R)) + 2nF 2 R + y) + \v\ 2 {F 2 - f 2 ) 2 R 3 

, (13) 

The additional complexity is that we must now minimize 
Eb for the single soliton energy with respect to F for each 
R or at each density. 

Further relaxing of constraints or solving the full set 
of Euler-Lagrange equations for each field is somewhat 
laborious, so this is as far as we go. This is so because 
this Wigner-Seitz approximation has other uncertainties. 
We refer the reader to a more sophisticated treatment of 
this problem pfij] . 

This section, as stated earlier, is to be seen as peda- 
gogical. It does not yield a realistic equation of state. 

Clearly the energy eigenvalue goes up as R (width of 
the potential) is decreased, till, at some lowest value of 
R, a bound state solution cannot be supported and this 
phase is lost. 

Remarks 

i) The Equation of State (EOS) with the above vari- 
ational ansatz is obviously inadequate to produce 
the single soliton/nucleon with arbitrary variations 
of < a > and < ir >. The exact solution fits the 
nucleon for a value of g y = 5.4. This will not be 
the case here. We will need a larger coupling, g y , 
to fit the nucleon in this approximation. Clearly, 
the minimum of E^ vs R, at i? m i n , is the energy of 
a single isolated soliton. 

ii) The nucleon may be obtained by projection of the 
soliton into good spin isospin states, J — I = ^. 
The soliton is largely a wave packet, a linear su- 
perposition of all J = I = (n + states. The 
maximum weight comes from the lowest, J = I 
states. We make the approximation that it is an 
equal linear superposition of the Nucleon (N) and 
A states and set the soliton energy to be midway 
between Mn and Ma- 

E b , min = M so iito„ = M N + - (M A - Mn) (14) 
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iii) Even with this approximation we find that to fit 
-E&,min> so as to produce a nucleon with Mn = 
940 MeV requires g y ~ 6.3. 

iv) Furthermore, there is a zero point quantum me- 
chanical energy of localization associated with the 
localized nucleon states that make up the crystal 
that may be estimated from the Uncertainty Prin- 
ciple. 

v) Another correction is that due to one gluon ex- 
change interaction between the quarks. 

We leave out the last two corrections and others as we 
shall be using a more realistic EOS anyhow (see below). 

Figure 1 shows the energy per baryon, Eb versus the 
baryon density rib in this phase. In this, the minimization 
of Eb with respect to the variational parameter F has 
been carried out. 

As we have pointed out, to construct the exact ground 
state of the crystal is, to say the least, a formidable ex- 
ercise. We have tried to employ an educated variational 
ansatz in the hope that it provides a good approximation. 
The EOS this yields has many correct qualitative features 
but nevertheless is much stiffer than most known EOS - 
it may not be unfair to say that it is almost arthiritic. In 
this instance it may be more judicious to use a well worn 
nuclear equation of state, like the APR98 EOS j^, for 
the entire nuclear phase, provided the density at which 
the transition to quark matter takes place is not much 
above twice nuclear density. A comparison between our 
stiff EOS and the APR98 is provided in Figure (1). 

The net result is that the nucleon is a soliton, the nu- 
cleons interact with each other to produce binding at 
around nuclear density. We then expect to have a Fermi 
liquid of interacting nucleons till the nucleons begin to 
overlap. As the density increases the Eb increases due 
to topological repulsion, we expect that the nucleons are 
no longer free to move around and get localized into a 
crystal (of nucleons) like configuration analysed above. 
This solitonic phase then dissolves into quark matter. 

IV. TWO FLAVOUR QUARK MATTER 

We shall now consider in Mean Field Theory the phases 
of two- flavour quark matter in the SU{2) l x SU(2)r chi- 
ral model above. 

We shall then extend the model to three flavours (u, 
d, s). 

A. The space uniform phase 

We now turn to the phase in which the pattern of sym- 
metry breaking is such that the expectation values of the 
meson fields are uniform. At zero density they are just 
the VEVs. 

<a> = U (15) 
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FIG. 1: Energy per baryon as a function of baryon density in 
the Wigner-Seitz quark soliton crystal as a model of nuclear 
state (solid line). The dashed line shows the corresponding 
relation in the APR98 pH equation of state for beta-stable 
matter, obtained using A18+<5«+UIX interaction model. 

< 7? > = (16) 

For arbitrary density we allow the expectation value to 
change in magnitude, as it becomes a variational param- 
eter that is determined by energy minimization at each 
density. 

< cr > = F (17) 

< 7? > = (18) 

Such a pattern of symmetry breaking simply provides a 
constituent mass to the quark m = g y < a > = g y F and 
the quarks are in plane wave states as opposed to the 
bound states in the nucleonic phase [l3[ . 

The mean field description of this phase is simple. The 
energy density is 

^ = E 7^37 / d^kVr^T^ +^«* 2 > -flf 
u,d ^ ' 

where m = g y < a >= g y F and the degeneracy 7 = 6. 
We shall use g y = 5.4 as determined from fixing the nu- 
cleon mass in this model at 938 MeV The integral 
above runs up to the 'u' and 'd' fermi momenta. 

For neutron matter (without /3 equilibrium) we have 
the relations 

H = (ir 2 n u )i = (n 2 n b )i (20) 
k{ = (2TT 2 n b ) L * (21) 

E b = ^ (22) 
rib 

where rib is the baryon density. 

At any density the ground state follows from minimis- 
ing free energy, with respect to < a >= F. As shown in 
the figures of 0,0,0], this phase begins at n& = 0, with 
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Eb = Sgfir, which then falls till chiral restoration occurs 
at some nx- After this, as the density is increased, 
continues to drop, goes to a minimum and then starts 
rising corresponding to a massless quark fermi gas. 

In the chirally restored phase the EOS is very simple 
and parallels the MIT bag description of 22] : 



n b > n x 



3 \ 2 I , A 2 „ 4 



(23) 
(24) 



The last term above is just the bag energy density, and 



a= (1 + 23) 



(25) 



This phase has two features, a) chiral restoration at nx 
followed, with increasing density, by b) an absolute min- 
imum in at a np > nx- 

From the comparison of this phase with the nucleon 
and nucleonic 'phase' arising from the same model (see 
El El) > it is clear that the nucleonic phase is always 
of lower energy than the uniform phase above, upto a 
density of roughly 3 times nuclear density, which is above 
the chiral restoration density in the uniform phase. 



B. The Pion Condensed phase 

Here we shall consider another realization of the expec- 
tation value of < a > and < 7? > corresponding to pion 
condensation. This phenomenon was first considered in 
the context of nuclear matter. 

Such a phenomenon also occurs with our quark based 
chiral a model and was first considered at the Mean Field 
Level by Kutschera and Broniowski in an important pa- 
per 0- Working in the chiral limit they found the pion 
condensed state has lower energy than the uniform, sym- 
metry breaking state (phase 2) we have just considered, 
at densities of interest. This is expected as the ansatz for 
the PC phase is more general than for phase 2. 

The expectation values now carry a particular space 
dependence 



< a > 

< 7T 3 > 

< 7Ti > 

< 7T 2 > 



F cos (q.r) 
F sin (q.r) 





(26) 
(27) 
(28) 
(29) 



Note that when |<f| goes to zero, we recover the uniform 
phase (2). The Dirac Equation in this background is 
solved in fijl and reduces to 



Hx(k) = (d.k- ^q.aj 5 T 3 



(3m) X (k) = E(k) X (k) (30) 



where m = g y F . The interaction term has been recast in 
terms of the relativistic spin operator, cr/s. It is evident 
that if spin is parallel to q and T3 = +1 (up quark) then 
this term is negative and if T3 = — 1 (down quark) then it 



is positive. For spin antiparallel to q the signs for T3 = +1 
and —1 are reversed. 

The spectrum for the hamiltonian is the quasi particle 
spectrum and can be found to be 



E(_)(k) = \/m 2 + k 2 + iq 2 - \Jm 2 q 2 + (q.k) 2 (31) 



E (+) (k) = ^m 2 + k 2 + -q 2 + ^m 2 q 2 + (q.k) 2 (32) 

The lower energy eigenvalue £<(-) has spin along q for 
T3 = 1 , or has spin opposite to q for T3 = — 1 . The higher 
energy eigenvalue has spin along q and T3 = —1, or 
has spin opposite to q and t 3 = +1. 

In this background the fermi sea is spin polarized into 
the states above. The quasi particles are, however, good 
states of T3. 

First we fill up all the lower energy, E(_^(k), states 
and then we have a gap and start filling up the Er + \ (k) 
states till we get to E F , the fermi energy corresponding 
to a given density for each flavour. 



1 



f7 



d 3 kO(E l F -£?(_)(*)) 



+ / d 3 kQ(E F -E (+) (k)) 



(33) 



rib 



(n u + n d )/3 (34) 
-^7 (J d?kE { ^(k)Q(E F - E^(k)) 



(2 



d 3 kE (+) (k)e(E F ~E {+) (k)) 



y(^ 2 -/. 2 ) 2 



(35) 



(36) 



We can now write down the equation of state as in 
It is found that the PC state is lower in energy than the 
uniform phase 2 for densities of interest. For the explicit 
numbers and figures we refer the reader to 0- 
We briefly remark on some features of this phase. 

1. The reason that the PC phase has energy lower 
than the uniform < a > condensate is perhaps 
best understood in the language of quarks and anti 
quarks. To make a condenste a quark and anti- 
quark must make a bound state and condense. For 
a uniform < a > condensate the q and q must have 
equal and opposite momentum. Therefore, as the 
quark density goes up the system can only couple 
a quark with k > kf and a q with the opposite mo- 
mentum. This costs much energy so the condensate 
can only occur if kf is small, at low density. On the 
other hand, the pion condensed state is not uni- 
form. So at finite density, if we take a quark with 
k = kf the q can have momentum k — \kf — q\, 
which is a much smaller energy cost 

2. Since the pion condensate is a chirally broken 
phase, the chiral restoration shifts from very low 
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density in the uniform phase to very high density, 
~ 10p lmc . This is a signature of this phase. 

3. Since this phase is lower in energy than the uniform 
phase for all densities of interest we go directly from 
the nucleonic phase to the PC phase completely by- 
passing the uniform phase showing that all the 'in- 
teresting' features and conjectures for the uniform 
phase are never realized. 

4. Another feature of this ttq condensate is that since 
we have a spin isospin polarization we can get a 
net magnetic moment in the ground state, as the 



magnetic moments of the u and d quarks add. 



V. THE THREE FLAVOUR STATE 

The extension to three flavours or SU (3) chiral sym- 
metry needs some clarification. 

The generalized Dirac Equation for the SU (3) case is 
considerably more complicated and involves a singlet £o 
and an SU(3) octet £ a of scalar fields and a singlet 4>o and 
an SU(3) octet <fi a of pseudoscalar fields, that interact 
with the quarks as shown in |23j. 



Hi>{k) = (-ia.d - 5 ^(v/2/3(Co + ^o7s) + A a (£ a + i<t>als)))i> = ^ 

I 



(37) 



In the chiral limit, the spontaneous symmetry breaking 
pattern is not unique. We choose the pattern in which 
the SU(3)l x SU(3)r chiral symmetry breaks down to a 
vector SU(3). For the uniform case, we have 



< Co > 

< > 

< 00 > 

< 0a > 



(38) 
(39) 
(40) 
(41) 



This gives a constituent mass m — g y f n for all (u, d and 
s) quarks. The explicit symmetry breaking strange quark 
mass term with current mass m s , is then added to H. 
The strange quark mass, M s , then, turns out to be the 
sum of the constituent and explicit mass, M s = g y fn+m s 

Here we run into a puzzle of sorts. The mass term now 
has two components, an explicit or current mass, m s and 
a constituent mass, g v < a >= g y j^. However, it does 
not disturb any relation (e.g. the GT relation), whether 
we choose, < a > = + or — f w . This raises an ambiguity 
about the relative sign of the two mass terms. It would 
seem that the choice of opposite signs for the two terms 
is optimal since it gives the lowest mass or ground state 
energy for the one fermion sector. 

When the strange quark gets a large explicit mass 
m s about 150 MeV (the u and d quarks have negligi- 
ble explicit masses), this question becomes very relevant. 
However, experimentally the strange baryons have larger 
masses than the non strange ones suggesting that the 
relative sign is plus. Georgi, in his book [24|, finds that 
the success of the non relativistic quark model, partic- 
ularly for the masses and the magnetic moments of the 
baryons, follows from taking constituent masses of about 
350 MeV for the u and d quarks and a total mass of 
about 550 MeV for the strange quark. Clearly, the larger 
mass of strange baryons suggests the same. Thus, exper- 
imentally, the relative sign, plus, is selected. We shall 
therefore continue to use this. 



A. The three-flavour Pion Condensed phase 

For describing strange quark matter we use the 3- 
flavour Pion Condensed state This is a more ver- 
satile state than the one used in [2^ (3 flavour Chirally 
Restored Quark Matter - CRQM) , the latter being a sub- 
set of the former. 

Next, we formulate the symmetry breaking in the pres- 
ence of the pion condensate. This is given as follows, 

< Co > = \/3/2F(l + 2 cos (<f.r))/3 (42) 

< > = -V3F(1 - cos (<f.r))/3 (43) 

< 0o > = (44) 

< 3 > = F{sm{q.r)) (45) 

while all other fields have zero expectation value. 

This gives exactly the PC hamiltonian equation for the 
u, d sector and continues to give the simple mass relation 
for the strange quark, M s = g y F + m s . When q = and 
m s = we recover the chiral limit above. 

We may now simply add the two-flavour PC results 
for the energy density and density derived above to the 
strange quark energy density which arises from the single 
particle relation, 



E, 



y/Mj + fc 2 



The strange quark energy density is given by Baym (eqn. 
8.20) H 



—Mj(x s n s (2x] 

07T 



+ 1) - ]n(x s + n s )) 



(46) 



where x s = k{ /M s and n s = *J\ + x 2 , k( being the fermi 
momentum of the strange quarks. 

The total energy density of the quarks for the 3 flavour 
PC is given by 



1 



A 



+ £ d + e s + -F\ 2 + ^-{F 2 -{fl)f (47) 
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From the effective potential given in 

[E for the 5/7(3) 
case, there is an extra factor of 3/2 that multiplies the 
last term. This can be absorbed, a s we have done, by a 
redefinition, Ai = AX, where A = a/3/2. 

The correction due to one gluon exchange interaction 
can also be incorporated in the evaluation of e u , e d and 
e s above. We have discussed this in some detail in 
In this paper we follow the prescription of Baym pEj , as 
done in to calculate the interaction contribution to 
quark energy densities. 



B. /3-equilibrium in the PC phase 

We have the following general chemical potential rela- 
tions for quark matter 



E F 
E d F 



/<>, 

Md = Us 



Me — fJ-d — Hu 

n - A 
e ~ 3tt 2 



(48) 
(49) 
(50) 

(51) 



The charge neutrality condition below further reduces the 
number of independent chemical potentials to one. 



2n u (n u ,q, F) - n d (pd, g, F) - n s (fi s ) 
3 

The baryon density is 



n e = (52) 



rib 



JjXu, q, F) + n d (n d , q, F) + n s (p s ) 



n s = (k{f/(^) 



(53) 
(54) 



For matter in [3 equilibrium we need to add the electron 
energy density to the quark energy density above 

e e = (1/4tt 2 )^ 
The total energy density is 



TABLE I: Charge neutral, 3-flavour, beta-equilibrium pion 
condensed phase with m CT = 800 MeV. The columns are: u- 
quark chemical potential (p u in MeV), baryon density (n& in 
fm -3 ), energy per baryon (Eb in MeV), electron density (n e 
in fm -3 ), ratio of densities of d-quark and u-quark (nj/n u ), 
that of s-quark and u-quark (n s /n u ), the order parameter (F 
in MeV) and magnitude of the vector q. 



flu 



lib 



E b 



F 



280.0 
300.0 
320.0 
340.0 
360.0 
380.0 
400.0 
420.0 
440.0 
460.0 
480.0 
500.0 



0.2972 
0.3645 
0.4591 
0.5409 
0.6700 
0.7628 
0.8967 
1.0472 
1.1774 
1.3225 
1.5149 
1.6900 



984.94 
981.48 
994.07 
1008.88 
1043.21 
1062.14 
1104.94 
1134.64 
1168.76 
1216.52 
1246.21 
1290.81 



.2303E-02 
.1602E-02 
.1141E-02 
.1173E-02 
.9455E-03 
.1063E-02 
.4413E-03 
.1283E-02 
.4890E-03 
.1794E-03 
.3896E-03 
.2212E-03 



1.916 
1.731 
1.599 
1.564 
1.499 
1.482 
1.345 
1.463 
1.317 
1.218 
1.267 
1.211 



.0318 
.2269 
.3640 
.4004 
.4672 
.4847 
.6245 
.5040 
.6529 
.7529 
.7033 
.7597 



37.0406 
31.9655 
28.6471 
30.4335 
28.8195 
31.6743 
23.6066 
36.1496 
26.8871 
19.3370 
26.9453 
23.0806 



2.5945 
2.6149 
2.9703 
3.0216 
3.6217 
3.4574 
3.9497 
3.9839 
4.0720 
4.3125 
4.3670 
4.4865 



The baryon chemical potential is defined as 

fi b = de/dn b (56) 

After meeting all the neutrality and equilibrium condi- 
tions above for fixed F and q, we can write all the above 
variables as a function of a single variable, fi u . We then 
minimize Q independently with respect to F and q. The 
Eb etc then follow. 

The results are presented in the tables below and in 
Figure 2a. Fig 2b gives the the EOS for all the phases 
considered so far. 



VI. THE PHASE DIAGRAM AND THE 
EQUATION OF STATE 



e = e p + e e 

The energy per baryon, Eb — e/nj,, then follows. 

For the pion condensed state, the ground state energy 
and the baryon density depend on the variational pa- 
rameters, the order parameter or the expectation value, 
F = V< 7? > 2 + < er > 2 and the condensate momentum, 
\q\. To define the free energy at a fixed baryon density 
then requires some care. 

We go about this by defining a baryon chemical po- 
tential to go with a baryon density. We have obtained 
both the baryon density and the energy density of the 
PC in terms of the u,d,s quark fermi energies/chemical 
potentials. First we construct the free energy 

£1 = e - ribfJ-b = e P + e e - nbUb (55) 



The starting point for the phase diagram of QCD at fi- 
nite density is as such: At very low density we know that 
there is chiral SSB, with the pion as the Goldstone bo- 
son - this breaks chiral symmetry spontaneously, leaving 
colour symmetry unbroken. At very high density we have 
a colour SC pairing instability for the quarks and of the 
many pairings investigated |26| the CFL (colour flavour 
locked) pairing is favoured - this spontaneously breaks 
colour and chiral symmetry. In between these limits the 
issue of the ground state is open. 

We approach this problem, from the low density end, 
by the effective lagrangian L that can access both, nu- 
clear and quark matter states, with or without chiral 
SSB. We have argued that this L has validity upto en- 
ergy scales of about 800 MeV. The advantage of this L is 
that it has only 3 couplings and no other free parameters. 



10 



Pion Condensed Matter with Strange Quarks 
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FIG. 2: (a) Upper panel: energy per baryon vs baryon num- 
ber density for 3-flavour pion-condensed phase for three val- 
ues of assumed tree-level mass of the scalar meson a. Charge 
neutrality and beta equilibrium are imposed. One-gluon ex- 
change interaction is included using the prescription of Baym 
|25| . (b) Lower panel: comparison of the equations of state of 
the 3-flavour space uniform phase and pion-condensed phase 
for m a = 800 MeV, with APR98 plj. 



These can all be determined from experiment. However, 
at the mean field level, this L is good for investigating 
chiral condensates but not colour condensates. 

If we start at the high density end, as done by the 
several authors who practice diquark colour supercon- 
ductivity, one knows the ground state at very high den- 
sity, when asymptotically free QCD makes calculation 
reliable. However, for intermediate density phenomeno- 
logical four fermion interactions have to be introduced. 
This is equivalent to introducing a given value for the 
gap parameter, A |26(. Also, another phenomenological 
input is the bag parameter, B. Uncertainty attends these 
parameters. It is not clear how this ground state transits 
into the low density chiral SSB state. 

We start with saturation nuclear matter at nuclear 
density and this persists till the nuclear matter gets 
squeezed into quark matter at moderately higher den- 
sity, when a pion condensed state takes over p|. In 
all these phases, as we have found, chiral symmetry is 



spontaneously broken - though the patterns of symme- 
try breaking keep changing with baryon density. Since 
we have not investigated all other possible condensates, 
we cannot vouch for our neutral pion condensate being 
the best ground state. 

It is. however, to be noted from the lecture notes of 
Baym |25| . that the neutral pion condensate is preferred 
over the charged pion condensate for charge neutral nu- 
clear matter, in the non relativistic limit (particularly if 
we put the axial coupling constant gA — 1)- We have 
found that this is the case for charge neutral quark mat- 
ter as well 27]. It is worth pointing out that all these 
states that have lower energy than the chirally restored 
CRQM state, are chiral symmetry broken states. 

At even higher density the most likely state is a diquark 
condensate - a colour flavour locked (CFL) |2^| state - 
which can persist till arbitrarily high density. Such a 
ground state spontaneously breaks both chiral and colour 
symmetry. Diquark condensates are unlikely at moder- 
ately high density as they depend on the quark density of 
states and so the pion condensate is most likely at such 
densities. 

There is an important issue that arises here and that 
is the comparison between the diquark condensate state 
and the pion condensed state. In this case the starting 
point is a chiral symmetric four fermion interaction which 
can accommodate both chiral (quark-antiquark colour 
singlet) condensates and diquark condensates. Such a 
comparison has already been done by Sadzikowski [2^ 
in the context of a NJL chiral symmetric model, for the 
case of 2 flavours - SU(2) L x SU(2) R . What is done is 
at the level of mean field theory. The NJL model has 
four fermion interactions in terms of the quark bilinears 
corresponding to the a and 7r field quantum numbers, 
with a common dimensional coupling, G. If we are in- 
terested in a ground state carrying sigma and/or pion 
condensates we can replace these quark bilinears by the 
corresponding a and tt EV's in the MFT. This yields 
the ground state energy of the space uniform SSB and 
the PC states. Alternatively, this NJL can be mapped 
to our linear sigma model and the ground state thereof, 
which has been considered earlier. 

To get to the diquark condensate state we have to Fierz 
transform this NJL chiral, L, and look for its projection 
into the diquark condensate channel and follow the same 
procedure of MFT. This projection gives the diquark con- 
densate lagrangian used by [2rl |2S| , with a four fermion 
coupling, G' = G/4@. 

Working with simultaneous MF condensates, corre- 
sponding to space uniform chiral SSB (which generates a 
spontaneous mass for the quarks (no PC)) and diquark 
SC (which gives rise to colour SC), they |2g find that for 
the above value of the two G s, that follow from the Fierz 
projection, the chiral condensate is always the preferred 
state up till the limit of validity of this model. 

However, they find that with arbitrary (higher) val- 
ues of G , it is possible to have a phase transition from a 
space uniform chiral SSB to a diquark condensate at some 
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large density. In particular, they consider, G = G, and 
find that from a pure chiral SSB state, with m = g y f v 
= 301 MeV, there is a continuous phase transition to a 
mixed chiral SSB with a small admixture of diquark con- 
densate state at a quark chemical potential of about 0.33 
GeV. For the higher value of, m — 500 MeV we employ, 
this is expected to occur at a higher value of fi. This is 
followed by a first order phase transition at slightly higher 
/i, to also a dominantly diquark condensate mixed state. 

This state then evolves to an almost entirely diquark 
condensate. 

Till now we have not cosidered the PC state. The 
PC state always has lower energy (free energy) than the 
uniform chiral SSB state. It would then be reasonable to 
expect that the PC state and not the uniform state would 
be the preferred state of chiral SSB. Since chiral SSB 
persists till much higher density (/i) in the PC state we 
expect that with the inclusion of this state the transition 
to the colour SC state will be pushed to higher density 

A following work by the same author addresses this 
matter by considering simultaneous MF condensates of 
the, PC, and the diquark condensate j2(|. It is found 
that the first order transition does shift to higher density 
(/j = 400 MeV). This is considered only for the partic- 
ular case G = G/2, and for m = 301 MeV. For our 
larger value of m this is expected to occur at larger \i. 
Furthermore, in this case it seems after the first order 
transition occurs, the chiral SSB order parameter, M, 
goes to zero, indicating that this state has no chiral SSB 
but only colour SSB. Of course, this is so as these works 
deal with the two-flavour case - where the colour diquark 
condensate is a chiral singlet. 

Realistically, we must consider 3 flavours, since the 
quark chemical potential is much greater than the strange 
quark mass. In this case we are very likely to have the 
CFL state as the lowest energy state. The criterion for 
this is given in [3(j and is A > m^/4/i, which is easily 
satisfied. Furthermore, in this case the diquark conden- 
sate is a colour-flavour condensate which has both chiral 
SSB and colour SSB, albeit in a manner different to the 
PC. 

The deciding question is then what effective L is to be 
used. There is no derivation of the exact effective L from 
QCD at arbitrary scales. There are many four fermion 
type (or higher order) interactions derived or rather mo- 
tivated from various points of view, e.g. instantons, gluon 
exchange etc. However, there is no unique L that is de- 
rived as such. It is therefore necessaryto have an ar- 
gument to decide this issue. See also [34| for a related 
discussion. 

Since the NJL chiral L identifies with our linear sigma 
model, and this latter model is what we have used as a 
valid model till centre of mass energies/scales of less than 
800 MeV, the right procedure would be to take this model 
to describe physics upto this scale. In this case, as we 
have argued, the PC is the preferred ground state, till the 
scale of validity of our effective L (this is for G = G/4). 



Even if we relax this to larger G we have a three-flavour 
PC, till n well above 400 MeV, followed by a CFL state, 
with increasing density. As the tables suggest such values 
of, /i, correspond to baryon density 5-6 times nuclear 
density. This makes the PC as likely state in neutron 
star cores. 

With increasing baryon density we then find the fol- 
lowing hierarchy. At nuclear density and above we have 
nuclear matter with chiral SSB, followed by the pion con- 
densed quark matter, again with chiral SSB, albeit with 
a different realization and finally a transition to the di- 
quark CFL state which also has chiral SSB (and colour 
SSB), with yet another realization. A point to note is 
that CRQM or free fermi seas are unstable. To one's 
surprise at zero temperature, at any finite density chi- 
ral symmetry is never restored! A similar result has also 
been obtained in 1+1 dimension 1351 . 



VII. DISCUSSION AND IMPLICATIONS FOR 
NEUTRON/STRANGE STARS 

A. Stars 

We have used the quark based linear sigma model for 
this analysis, where the tree level sigma mass was set 
to be around 800 MeV as indicated by the analysis of 
Schechter et al 14], by matching the results for meson- 
meson scattering from this model, with experiment. 

Such a tree level sigma mass also ruled against conven- 
tional SQM being the ground state of matter Strange 
stars require the absolute stability of SQM. Since this has 
been shown to be highly implausible so is then the exis- 
tence of strange stars. 

Independently, we find (Bhattacharya and Soni, in 
preparation) that realistic Neutron stars with PC cores 
occur only when the tree level sigma mass in this model 
is in a small window, 750 - 850 MeV. 

This can be seen from vs 1/nb diagram for the 
phases we have considered (Fig 3). The negative slope 
in this figure gives the pressure and the intercept on the 
vertical axis the baryon chemical potential. The Maxwell 
contruction of a common tangent to these curves gives 
the pressure and baryon chemical potential at which the 
transition occurs with a density discontinuity. 

We find that at m a = 800 MeV we get a transition 
with a relatively small pressure and a small density dis- 
continuity. This permits us to have a well developed PC 
core in a neutron star. The small density jump, besides, 
reassures us that a mixed [3lj phase would give results 
that would be very similar. Of course, with the APR98 
pl| EOS that we use, we do not have analytical expres- 
sions for all variables, which precludes a mixed phase 
calculation. 

At m a = 850 MeV, the transition moves to higher 
density and higher pressure with a larger density jump, 
so much so that there there is no PC core for the star - as 
the maximum mass instability for the star occurs before 
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FIG. 3: The Maxwell construction: Energy per baryon plot- 
ted against the reciprocal of the baryon number density for 
APR98 equation of state (dashed line) and the 3-flavour pion- 
condensed (PC) phase, for three different values of m CT (solid 
lines). A common tangent between the PC phase and the 
APR98 phase in this diagram gives the phase transition be- 
tween them. The slope of a tangent gives the negative of the 
pressure at that point, and its intercept gives the chemical po- 
tential. As this figure indicates, the transition pressure moves 
up with increasing m a , and at m a below ~ 750 MeV a com- 
mon tangent between these two phases cannot be obtained. 



a core can form. 

On the other side, at m„ = 750 MeV, the curves for the 
two phases do not permit a common tangent construction 
any more - so again a star with a PC core is ruled out. Is 
this a coincidence that a single parameter in our effective 
L, the mass of the sigma or A, plays a crucial role? Is it 
fortuitous that the tree level sigma mass set by scattering 
experiments sits in a small window that simultaneously 
rules out SQM as the absolute ground state of matter 
and also can provide us with neutron stars that can have 
pulsar range magnetic PC cores? 

The problem in sustaining a PC core with a nuclear 
exterior is that we have a stiff exterior with a soft interior 
- a rather unstable situation. It is thus not so surprising 
that very particular conditions must obtain for this to 
occur. 



B. The Effective Lagrangian and Mean Field 
Theory 

The question of the validity of our effective L is of 
the essence. Furthermore, since we do not go any further 
than MFT, the question of the validity of MFT is another 
equally important question. Linked to this is the question 
of QCD corrections. 

1. Meson sector 

Let us begin with the part of the Lagrangian that de- 
scribes the meson sector. As we have said earlier the 
justification for this, its range of validity and the value of 



the coupling A or the sigma mass have been very clearly 
spelt out in 01 . It will be useful to review this. As 
stated in this work it is the tree level Lagrangian that 
is used to describe 7r-7r scattering in the scalar channel. 
This tree-level amplitude is then improved by unitarising 
the K matrix. A very good fit to scattering data in this 
channel follows if we choose the sigma mass in the tree 
level L to be 800-850 MeV. The validity of this vis a vis 
the data is upto ^/s < 800 MeV. 

The actual mass and width of the sigma can then be 
gleaned by looking at the pole in this K unitarized am- 
plitude and is found to be 460-600 MeV. 

We have used this to fix the tree level coupling/sigma 
mass for our L. Since we work only in MFT (tree level), 
this also perhaps tells us that we may be off the mark by 
about 30%. 

Another cautionary point is that the above results are 
based only on scalar channel scattering and may not ap- 
ply generally. It is also possible to fit the data with chiral 
perturbation theory (infinite m a ) and the p resonance. 
It is possible that by excluding the p we may be miss- 
ing some short-distance repulsion, which might result in 
making our PC EOS too soft. 

2. Quark sector 

For obvious reasons one cannot find an analagous scat- 
tering experiment for the quark meson sector. However, 
there is other indirect evidence in this sector. 

i) The quark soliton nucleon that we have used to fix 
the coupling, g y , satisfies the Goldberger Tricmann 
relation exactly [3^J in accord with chiral symme- 
try. As stated in the introduction several indepen- 
dent properties of this nucleon, for example, mag- 
netic moments compare well with experiment. 

ii) In the large N c expansion, the soliton (and other 
physics) receives 1/N C corrections from loops. 

iii) Goksch ^2] finds that the finite temperature 
screening masses of mesons computed in lattice 
QCD compare very favorably with those calculated 
using finite temperature MFT for precisely our L 
(without including gluons). However, the couplings 
used are slightly different to ours, g y — 3.3 and 
sigma mass of 600 MeV. 

Georgi and Manohar Q use a non linear sigma model 
version of the above L to do effective cutoff field the- 
ory and argue for larger energy scale of chiral symmetry 
breaking vis a vis confinement. Using this scale for the 
cutoff, they are able to naively but consistently include 
arbitrary meson/quark loops. Further, they heuristically 
argue that the QCD coupling is weak for the chiral SSB 
vacuum and thus gluon loops may be ignored. 

Actually, as we have stated in the introduction, QCD 
can have multiple scales - a confinement scale, a chi- 
ral symmetry restoration scale and a compositeness scale 
for the pion. We would like to point out that the argu- 
ments of would go through if the compositeness scale 
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was substituted for the chiral restoration scale, since this 
scale is used for a momentum cutoff. This will allow the 
chiral restoration scale to fall anywhere in between the 
confinement and compositeness scale - a likely possibility 
suggested by the lattice. 

Parenthetically, we remark that the non-perturbative 
input of a finite mass sigma particle/resonance may do 
the same job. This is also supported by observations in 
Ref. and references therein. This also supports our 
effective L, which is built on the premise that the chi- 
ral SSB energy scale is larger than than the confinement 
scale. 

3. QCD sector 

We have no evident justification for neglecting this strong 
sector or at most doing simply one gluon exchange cor- 
rections. 

Ref. |3| makes the heuristic argument that the gluon 
coupling becomes strong to precipitate chiral SSB, after 
which it is expected to be perturbative in the new bro- 
ken vacuum. They also argue, that from the colour and 
electromagnetic hyperfine splittings in the baryon sec- 
tor of the non relativistic quark model one may infer an 
a s = 0.28. 

At finite density, there is screening in quark matter. 
The correct expansion parameter is then the screened 
charge and not the usual running coupling constant writ- 
ten as a function of the quark chemical potential which 
is commonly used. Such a screened charge has been for- 
mally constructed by one of us [3i| and is clearly much 
smaller than the usual running coupling. This makes a 
perturbative expansion in this coupling more plausible 
for dense quark matter. 

Further, the ground state of quark matter, the pion 



condensate, has spin-charge ordering. Long range forces 
like one gluon exchange then inhibit any change in the 
ground state and stabilize it against fluctuations as in 
many condensed matter situations where the coulomb 
forces are operative. In this case MFT may be valid. 

All the calculations referred to above for quark matter 
states have been carried out in MFT. We have tried to 
give some justification for its use. However, we cannot 
provide any rigorous proof for the validity of MFT. It 
is also the simplest thing to do and works well in many 
cases even if the coupling is strong. 

Since we began this work with an eye to neutron stars 
it may be appropriate to present our finding. We expect 
that the density profile of the star will start with the 
nucleonic EOS on the outside and go to a pion condensate 
in the interior and could well go to a colour flavor locked 
state at the centre if the density there is large enough; 
though, this is unlikely if we have a pion condensate as 
its softness may result in a smaller maximum mass. 
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